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Abstract 

We calculate the s-wave pion-pion scattering length in the isospin 1 = 2 channel in 
lattice QCD for pion masses ranging from 270 MeV to 485 MeV using two flavors of 
maximally twisted mass fermions at a lattice spacing of 0.086 fm. Additionally, we 
check for lattice artifacts with one calculation at a finer lattice spacing of 0.067 fm. 
We use chiral perturbation theory at next-to-leading order to extrapolate our results. 
At the physical pion mass, we flnd rriT^al^"^ = —0.04385 (28)(38) for the scattering 
length, where the flrst error is statistical and the second is our estimate of several 
systematic effects. 
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1 Introduction 



In the limit of massless up and down quarks, the resulting chiral symmetry of 
QCD is spontaneously broken and consequently the meson spectrum contains 
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three massless Goldstone bosons, vr^ and vr*^. Due to the unique role played by 
the pions, their interactions are strongly determined by the underlying chiral 
symmetry, and the s-wave pion-pion scattering lengths even vanish in the chiral 
limit. Thus the scattering lengths are sensitive to the chiral dynamics of the 
strong interactions, and the non-perturbative calculation thereof, the subject 
of this paper, is an integral part of understanding the low energy properties 
of QCD. 

In nature, the masses of the quarks are not zero but small and induce an 
explicit but weak breaking of chiral symmetry. Correspondingly, the pions 
are not massless but light. This breaking of chiral symmetry is systematically 
treated in chiral perturbation theory (xPT) by considering the quark masses as 
perturbations. Furthermore, the pion-pion scattering lengths no longer vanish 
and at leading order (LO) in xPT are predicted by Weinberg [1] solely in 
terms of the pion mass, rriT,-, and the pion decay constant, /,r, as 

mX=°~TTr^ = 0.160(1) and m^^f ^ --^ = -0.0456 (1) , 



where a^^" and a^^^ denote the isospin 7 = and 1 = 2 s-wave scattering 
lengths respectively. The next-to-leading order (NLO) corrections depend on 
unknown low energy constants, which can be determined from experimental 
measurements or lattice calculations. 

The experimental measurement of — 7r+7r~e^z/ (Ke^) decays by E865 at 
BNL [2] gives 

m^a[=° = 0.203 (33) and m^a^;^ = -0.055 (23) . 



When combined with constraints from xPT, these measurments yield 
m^aif = 0.216 (14) and m^aif = -0.0454 (34) . 

A combination of several experimental and theoretical inputs from CGL [3,4] 
produces a consistent but more precise result of 

m^aif = 0.220 (5) and m^aif = -0.0444 (10) . 

Additionally, the recent measurments of K^a decays [5] and — )■ ir^n^ir^ 

decays [6] by NA48/2 at CERN [7] give, without making any use of xPT 
constraints, 

™ ^^=0 n 001 ™ ^^=2 



0.221 (5) and m^ai=^ = -0.0429 (47). 
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Including xPT in their analysis, NA48/2 finds [8] 

m^4=° = 0.220 (3) and mW^f = -0.0444 (9). 

The results are all consistent with each other and the most precise results from 
NA48/2 are in agreement with the lattice results given shortly. 

The non-perturbative determination of the scattering lengths from lattice 
QCD is possible, despite the Euclidean nature of the calculation, by using 
a finite size method due to Liischer [9-13]. This method capitalizes on the 
relationship between the energy eigenvalues of a two pion system enclosed in 
a finite spatial box and the scattering phase of two pions in infinite volume. 
In the derivation of this result, it is assumed that the physical box size is 
large enough to avoid significantly altering the two pion interaction. In this 
limit, the intrinsic finite size effects of each individual pion are exponentially 
suppressed. The dominant contribution to the finite size dependence of a two 
pion state is then simply given by the interaction between the two pions due 
to the finite volume. Thus the strength and nature, repulsive or attractive, 
of the interaction of two pions will shift the energies of the otherwise free pi- 
ons. This shift then relates the energy eigenvalues to the scattering phase and 
ultimately the scattering lengths of two pions. 

One obstacle to the lattice determination of the pion-pion scattering lengths 
is the presence of disconnected diagrams that render the calculation of the 
1 = channel computationally demanding. On the other hand, the simpler 
1 = 2 channel does not require such diagrams and consequently many lattice 
groups have focused their efforts on this case. Furthermore, most calculations 
of the scattering lengths to date have been carried out within the quenched 
approximation [14-35]. There have been only two previous calculations of 
a^^^ with dynamical fermions. The first such calculation was performed by 
CP-PACS with Nf = 2 tadpole-improved clover fermions at rather heavy pion 
masses in the range = 0.5 GeV to 1.1 GeV [36]. However, it is doubtful 
that xPT at NLO, or any order, can be applied to such heavy pion masses. 
The other full QCD calculation was performed by NPLQCD with domain-wall 
valence quarks on the Nf = 2 + 1 asqt ad- improved coarse MILC ensembles 
with = 290 MeV to 590 MeV [37,38]. Mixed-action xPT at NLO was used 
to perform the chiral and continuum extrapolations. At the physical pion mass, 
NPLQCD finds 

mXf = -0.04330 (42) and = ^,phy) = 6.2 (1.2) , 

where /^^^(/x) is a low energy constant (LEC) appearing in the xPT description 
of the quark mass dependence of the scattering length. As discussed later, 
/^^^(/x) is evaluated at fj, = f-K,phy, where /,r,phy is the physical value of the 
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pion decay constant. 

In this work we determine the s-wave 1 = 2 pion-pion scattering length and the 
corresponding Z^^^ by using the Nf = 2 maximally twisted mass fermion en- 
sembles from the European Twisted Mass Collaboration (ETMC). Our lightest 
pion mass is lighter than those of the previous calculations and allows us to 
further probe the chiral limit. Due to the properties of twisted mass fermions 
at maximal twist, our calculation is automatically accurate to O(a^) in the 
lattice spacing, a. Additionally, we perform an explicit check for large lattice 
artifacts with a single calculation at a finer lattice spacing. As presented later, 
we find at the physical pion mass 

mXf = -0.04385 (28) (38) and /^f (/i = /,,phy) = 4.65 (.85)(1.07) , 



which is in agreement with the above experimental measurements and phe- 
nomenological analysis as well as the previous lattice calculation. 



2 Method 



2.1 Liischer's finite size method 



As mentioned in the introduction, Liischer's finite size method relates the 
energy levels of two pion states in a finite volume to the scattering phase in 
the infinite volume. For the case of two pions with zero total three-momentum, 
this method establishes a relationship between the lowest energy eigenvalue 
E^^ with a given isospin / in a finite box of size L and the corresponding 
scattering length a^^. For the 1 = 2 channel, it is given in Ref. [10] as 



SEif = Eif-2m^ 



1=2 



47ra^^ „ 



J=2 



1 + ci- 



L 



+ C2 



J=2' 



L 



(1) 



where Ci = —2.837297 and C2 = 6.375183 are numerical constants. Thus the 
above result allows us to convert a lattice determination of the energy shift, 
SE^^"^^, into a calculation of a^^^. 
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2.2 Extraction ofSE^-'^ 



To extract SE^^"^, we construct the 7r+ and tt+tt^ two-point correlation func- 
tions from the operators proposed in Ref. [24], 

C^{t) = {{n+y{t + t,)n+{ts)) 

and 

a.(t) = ((vr%+)^(t + t.)(7r+7r+)(t,)). 

Here tg is an arbitrary time shce, vr~'"(t) = J2xid'~f5u){x,t) is an interpolating 
operator for the vr"*" meson with zero total three-momentum and (7r+7r'^)(t) is 
an interpolating operator for the two pion state, again with zero total three- 
momentum, given by 

(7r+7r+)(t) = 7r+(t + a)7r+(t) . 

In order to avoid complications due to Fierz rearrangement of quark lines as 
discussed in Ref. [24], we use the 7r+ interpolating fields at time slices separated 
by one lattice spacing. 

From the large time behavior of C^(t) and C7r7r(i), it is possible to extract the 
corresponding ground state energies as follows, 

C^{t) exp(-m^ t) and C^^(t) A^^ exp(-E^=^ t) , 

where we assume that t is large enough to neglect excited states but still far 
enough from the boundaries to ignore boundary effects. Furthermore, con- 
structing the following ratio of correlation functions we can determine 
directly as 

CTTTri^) ^TTTT / r- 7-t/=2 j.\ 

— ^—eM-6E„ t) 

where t satisfies the same requirements as before. However, we use anti- 
periodic boundary conditions for the quarks in the time direction in order 
to match the sea quarks used in our calculation, and this leads to a more 
complicated time dependence for and Cj^^. 
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2.3 Anti-periodic boundary conditions 

As mentioned above, in our calculation we employ anti-periodic boundary 
conditions in the time direction for the fermions. Using the transfer matrix 
formalism, the time dependence of our correlation functions is given by 

(Ot(t)O(O)) = Tr (e-^(^-*)O^(0)e-^*O(0)) / Z , 

where the time-slice transfer matrix is e""^^, the partition function Z is given 
hy Z = Tr(e~^"^) where T is the total time extent of our lattice and 0{t) 
represents either vr"'"(t) or (7r"'"7r"'")(t). Inserting a complete set of eigenstates 
of H into the above equation yields 

(Ot(t)O(O)) = ^ I (n|0|m) pe-^'-C^-^^e-^-V Z 

m,n 

= J2 I (ri|0|m) |2e-(^-+^")^/2 cosh((E^ - - T/2))/ Z . 

m,n 

The terms in the above series are thermally suppressed by factors of e"^""^ 
or e~^""^. Only those terms with Em = or = remain in the zero tem- 
perature, T —J- oo, limit. However, the effects of the suppressed contributions 
can still distort the behavior of correlation functions for finite values of T, 
particularly in the large t region. 

This phenomenon does indeed occur in this work for the two pion operator. In- 
termediate states {n\ = (vr+l and {m\ = {n~ \ give a constant, in t, contribution 

to C-JZTT) 

|(7r+|7r+7r+|7r-)|2e-''"-^/Z. 

This is comparable to the standard contribution, 

|(7^+7^+|7r+7^+|^])|2e-^-'^/2 cosh(E;^=2(t - T/2))/Z, 

when t approaches T/2. To be precise, for large enough volumes E^^"^"^ = 
2m^ -|- ~ 2m^, and hence these two contributions to Ctttt, e"'"'""^ and 

(.-^nTv T/2 cosh(£'7r7r(^ — r/2)) are in fact nearly equal for t = T/2. Additionally, 
the factor CT^{tY has similar problems. The correlator CT^{t) itself has a simple 
spectral representation. However, the square is more complicated and also 
contains a constant, in t, contribution as well. 
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To eliminate these contaminations, we use the derivative method [39] and 
define a modified ratio, R{t), in the following way 



The asymptotic form for R{t), ignoring terms suppressed relative to the leading 
contribution, is 

R{t + a/2) = An (cosh{6Eif t') + smh{6Eif t') coth(2m^t')) (3) 



where Ar is a combination of amplitudes in C-,^ and Cjr^ and t' = t + a/2 — T/2. 
Since m^r is the most accurately calculated component of our calculation, R{t) 
provides a nearly direct determination of SE^"^^ and cleanly eliminates the 
unwanted thermal contributions that spoil the simple ratio given earlier. 



3 Lattice Calculation 



3. 1 Twisted mass fermions 



In this work we use the two fiavor maximally twisted mass fermion configura- 
tions from ETMC Using twisted mass fermions at maximal twist ensures that 
physical observables are automatically accurate to O(a^) in the lattice spac- 
ing [40] . Most of the results presented here are from a sequence of ensembles 
with a lattice spacing of a = 0.086 fm and a box size of L = 2.1 fm. The pion 
masses range from 771.,^ = 270 MeV to 485 MeV. For the lower pion masses the 
volume is increased to L = 2.7 fm, and there is one calculation using a finer 
lattice spacing of a = 0.067 fm. The parameters relevant to this calculation 
are given in Tab. 1, and further details can be found in Refs. [41-43]. 



3.2 Stochastic sources 



For the calculation of pion correlation functions, it is known that the stochastic 
source method is more efficient than the point source method. Therefore, in the 
present work, we employ stochastic noise with two noise sources generated 
on each source time slice. Since we place the source on two time slices for the 
TT+TT"'" correlation function, tg and is + a, we therefore perform four inversions 
for each configuration. We remark that we also use the one-end trick in this 
work for the evaluation of correlation functions [44-46] leading to a further 
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Table 1 

Ensembles used in this work. Only dimensionless quantities are needed in this cal- 
culation, but for guidance we give the value of rounded to the nearest MeV for 
each ensemble indicated by /3, a/i and L/a. We also list the ratio mjr//^, the num- 
ber, A^, of configurations used, the energy shift a6El'^'^ and the scattering length 
myra^^^. The first uncertainty is statistical and, when present, the second one is 
systematic. 

improvement in the signal-to-noise ratio. Additionally, the source time slices, 
ts, are chosen randomly to reduce the autocorrelation between consecutive 
trajectories. 



4 Results 

4-1 Calculation of rriT^a^^'^ 

In Fig. 1 we show our lattice results for -R(t), defined in Eq. 2, as a function 
of the time t together with a correlated fit to the asymptotic form given in 
Eq. 3. All the ensembles shown in Fig. 1 visibly agree with the corresponding 
fit and lead to reasonable values of P^^ degree of freedom (dof), where 
is the correlated figure-of-merit function. To further verify these fits, we 
examined several possible sources of systematic error. First, the ratio could 
suffer from bias at large t, so we examined the jackknife estimate of bias 
but found it to be significantly smaller than the errors for all the ensembles. 
Second, we considered the possibility that we underestimated the errors due to 
autocorrelations. However, both the gamma method [47] and standard binning 
showed no significant signs of autocorrelation for R{t) in any of the ensembles. 
The possibility of vr*^ mixing, due to the breaking of parity at non-zero lattice 
spacing for twisted mass fermions, is considered in Sect. 4.2. But as described 
there in more detail, we find no statistically significant indications of the 7r° 
contributions. 
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71 




Fig. 1. The ratio R{t) as a function of t. The sohd lines are correlated fits to Eq. 3, 
from which the energy shifts aSE^"^"^ are extracted. The ensembles have been shifted 
vertically to facilitate easier comparison. 



There is one further possible systematic error due to the contributions from 
excited states in the small t region or from unphysical 7r° states in the large t 
region. To ensure that the fits for these ensembles are safe from such effects, 
we study the systematic errors caused by choosing a fitting window in which to 
match to the asymptotic form for R{t). First we ensure that the results exhibit 
clear plateaus when we increase the minimum t or decrease the maximum t 
used in the fits. However, to provide a quantitative estimate of the systematic 
error, we perform the following distribution method. We collect the results for 
aSE^"^^ from all fitting intervals with x^/dof < 2. This includes varying both 
the minimum and maximum time extent for the fitting range and results in 
30 to 60 values of aSE^^^ for each ensemble. We then make the distribution of 
these selected results and choose the median of this distribution for the central 
value. Then we take the central, and symmetric about the median, 68% region 
of the distribution to define the systematic error. Finally, we use the jackknife 
method to determine the statistical error on the central values. This method 
is also applied to m^a^^^, and the results for a5-E^^^ and m^a^^^ are given 
in Tab. 1. As shown in this table, the resulting estimates of the systematic 
errors are typically smaller than the corresponding statistical errors, and are 
at worst of the same order as the statistical errors. Since the distribution 
method used to estimate the systematic errors is itself subject to statistical 
errors, this is precisely what is expected if there are no substantial systematic 
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effects. However, since the final statistical precision for the value of mTrO^^^ at 
the physical limit turns out to be quite small, we decided, in order to avoid 
underestimating our final error, to carefully propagate these systematic errors 
through to the final result as described later in Sect. 4.5. 



4-2 71 contamination 

Twisted mass fermions violate parity and isospin at non-zero values of the 
lattice spacing. Therefore the spectral representation of the vr"*" and 7r'^7r+ 
correlators can admit states that would not be present in the continuum limit. 
In particular, unphysical contributions from the ir^, which has a mass m^o 
different, and smaller, than the mass of the vr^, may enter the Cj^ and C,r7r 
correlators in several ways [46]. Furthermore, these effects are believed to be 
more noticable in pion-pion scattering, so the successful calculation of all three 
isospin channels, / = 0, 1 and 2, would test the twisted mass formulation of 
lattice QCD. 

The 7r° can enter the correlator through intermediate states of the form 
(7r+|7r+|7r°) and (ir'^ir^lir'^lQ) . The former contribution is thermally suppressed 
by a factor of e~"^^°'^, however it leads to a time dependence with an energy 
of m^r — m^o that is lighter than the usually expected ground state. The 
second contribution is not thermally suppressed but corresponds to the first 
excited state with energy E^^+t^o ^ + m^o. This is lighter than the first 
physical excited state with energy near 3m^. Similarly, Ctttt contains unphys- 
ical contributions from (7r+7r+|7r"*'7r'''|7r'') and {7i~^'K~^7T^\7i^7r~^\Q) . Again there 
is an additional light state that is thermally suppressed by 6'"^^°"^ but has an 
energy of Et^+t^+ — m^o ^ 2m^ — m^o that is lower than the physical ground 
state near 2m^, and the first excited state is lowered to £'^+^+^o 2m7r + m^o 
rather than the expected energy of approximately 'i\l'm\ + [271 / LY. 

The parity violating matrix elements responsible for these effects are 0{a) in 
the lattice spacing, even at maximal twist, however the matrix elements appear 
squared in the correlators. Therefore these unphysical states make an O(a^) 
contribution. The question, however, is not about the scaling in the lattice 
spacing, but about the size of this contribution at the lattice spacings used in 
this work. A detailed discussion of this issue can be found in Ref. [48]. Here, 
our focus is more practical. We want to ensure that the scattering lengths 
calculated in this work are not significantly distorted due to these effects. 

First, the naive estimate for the suppression factor for the additional light 
contributions, — 771,^0 in C,r and Et,+^+ — m^o in Cnn, is (aAQCD)^e~'"'^°^. 
The value of m^o is difficult to calculate precisely, but it is clear from Ref. [48] 
that m^o is never more than 20% lighter than for the ensembles in this 
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work. Therefore we will simply use m^r and a value of Aqcd = 250 MeV to 
set the order of magnitude for these suppression factors. We find that for the 
ensembles used here, the largest value of (aAQCD)^e~"^TO^ is 9 ■ 10~^ for the 
(3 = 4.05, afi = 0.0030 ensemble in Tab. 1. Using the actual value of m^o 
from [48] raises this to 2 ■ 10~^. This value is small, but it is not too far beyond 
the statistical precision of the correlators used to calculate adE^"^"^, hence we 
must carefully check for these contributions. 

Second, there are the additional states that are only suppressed by (aAqco)^- 
However, these states are heavier than the physical state and hence would 
occur in the correlators as excited states. The naive suppression factors are 
1 ■ 10~^ and 7 ■ 10~^ for a = 0.086 fm and 0.067 fm respectively. These simple 
estimates are larger than for the other states, however these contributions are 
also more strongly suppressed by their own energies. 

In the light of these arguments, we made a significant effort to attempt to 
find such effects anyway. We tried fitting the individual C^(i:) and CT^nit) 
correlators as well as the ratio R{t) to various functional forms including the 
physical state and both the additional heavier and lighter states, just the 
lighter state or just the heavier state. We fit the most general forms, keeping 
all energies as free parameters, and additionally constrained forms, in which 
we constrained m^o based on known values. And we also explored several 
minimization methods. The net result was that one could indeed lower the 
value for each fit, but the P^r degree of freedom still increased, indicating 
no statistically significant contribution from the unwanted 7r° states. 

However, we must offer a few words of caution. While we could not find any 
compelling evidence for these contributions, we of course can not rule out their 
presence at a level beneath our statistical resolution. We should further note 
that there are visible excited states in the correlators. However, the accuracy of 
the correlators for the ensembles studied here does not allow us to distinguish 
the physical excited states, near 3m^ for C,r and 2^ml + (271 /LY for Ctttt, 
from the unphysical excited states, near + m^o 2m.^ for and 2m^ + 
m^o ^ SuItt for Ctttt- The extensive study of systematic errors due to the fitting 
range discussed in the previous section was partially motivated by these issues. 
It provides the quantitative statement that these effects do not rise to the level 
of our statistical precision and gives an estimate of the systematic error. 

Additionally, there are two reasons that these contributions may be smaller 
than anticipated. First, the unphysical contributions correspond to scattering 
states that may be suppressed by a power of the volume. Second, the construc- 
tion of R{t) in Eq. 2 forms a discrete approximation to the ratio of derivatives 
of C^^ and and may further suppress the nearly constant light state con- 
tributions. Finally, these effects are anticipated to be more substantial for the 
other isospin channels, hence a detailed understanding of the 7r° contribu- 
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tions to pion-pion scattering will have to await our ongoing calculations in the 
1 = 1 [49] channel and our planned work for 1 = 0. 



4-3 Finite volume effects 



The dominant finite size effect in this calculation is, of course, the shift in 
SE^"^"^ due to the interactions of two pions in a finite volume. Additionally, 
there are the exponentially small, as opposed to the merely power suppressed, 
finite volume corrections to / = 2 pion-pion scattering that have been de- 
termined for scattering near threshold in Ref. [50]. The resulting finite size 
corrections for the scattering length are given there as. 



where 



ml i ml d , ^, , 2m'i . , 
^Fv = -^^^AI{m^) + -f^AJ,.^,{Aml) 

J TT \ J TV TV J TV 



m^ 



213/2^5/2 1^ ) E ^rV—-lY 8 |n|m.L 



Oil 



^4 ^ g-|n|m,L ^ ^7 
|n|5^0 



Using the above result, we calculate the corrections to m^a^^^. Compared 
to the statistical errors, the finite volume corrections are negligible. To be 
precise, they are never more than 6% of the corresponding statistical error 
and are hence ignored in the following analysis. 

There is a second finite size effect originating from the effective range approx- 
imantion, ptan~^5(p) = + |reffp^, which is used to relate the scat- 

tering phase 5{p) at vanishingly small momentum p to the scattering length. 
The dependence on the effective range r^.^ is very small and gives rise to the 
corrections at 0{L~^) in Eq. 1. As argued in Ref. [38], assuming that the 
effective range is at most twice the scattering length, this correction can be 
estimated using the measured values of and SEl"^"^. Using the result given 
in Ref. [38] , we calculate this correction and find that it is never more than 9% 
of the corresponding statistical error of m^a^^^. Hence, this finite size effect 
is also sufficiently small to be ignored as well. 
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4-4 Lattice artifacts 



Most of the calculations presented here use a single lattice spacing of 0.086 fm, 
but we have also performed an additional calculation of SE^^^ and m^a^^^ at 
a second lattice spacing of 0.067 fm and at a pion mass of 307 MeV. This pion 
mass lies very close to that of the a = 0.086 fm, = 309 MeV point. The 
physical volumes of these two ensembles differ, so the values of can not 

be directly compared. However, assuming that Liischer's method correctly ac- 
counts for the finite volume dependence of SE^"^^ for these two ensembles, we 
can compare mTrO^^^ for the two lattice spacings, and indeed we do find statisti- 
cal agreement between the two ensembles as indicated in Tab. 1. Furthermore, 
as described in the next section, we note that the expected O(a^) corrections 
from maximally twisted mass lattice QCD are actually weakened to O(m^a^) 
for the I = 2, I3 = ±2 channel as shown using twisted mass xPT [51], thus 
suggesting further that the lattice spacing dependence of 772,^0^^^ is mild for 
the calculations in this work. 



4.5 Chiral extrapolation 



The pion-pion scattering lengths have recently been calculated in twisted mass 
xPT [51]. This is an expansion of twisted mass lattice QCD in both the quark 
masses and the lattice spacing. There it is shown that at NLO the lattice 
spacing corrections to the J = 2, J3 = ±2 scattering lengths are proportional 
to cos(co'), where uj is the twist angle. Thus at maximal twist, u = n/2, the 
explicit discretization errors vanish exactly, and the scattering length can be 
simply represented by the continuum NLO xPT formula [52,53]. 



,/=2 



As suggested in Refs. [37,38], we perform the chiral extrapolation of m,, 
in terms of mT^/ instead of m^r. Additionally, the xPT renormalization scale 
is fixed as /i = f-K,p\iy The resulting NLO expression is then 



1=2 



mz 



mz 



J TT 



1 - Ctt^IA* - /TT.phy) 



(4) 



where /^^ (/i) is related to the Gasser-Leutwyler coefficients k as [54] 

lifi^^) = 1^1 + ^4^2 -k- 4/4 + 3 In . 

3 3 jj,'^ 



It is important to note that extrapolating in m^^/ f^^ instead of simply 
does indeed change the expression for m-^al^"^ but only at the next-to-next-to- 
leading order (NNLO). The advantage of this form is that m-^j f.„ is calculated 
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directly on the lattice with small errors and the chiral extrapolation does not 
require fixing a physical value for the lattice spacing. 



We now fit our lattice results for m^a^^^ from Tab. 1 to the functional form in 
Eq. 4 in order to extrapolate m^a^^^ to the physical point and also extract the 
low energy constant /^^^(/i = /7r,phy)- The calculated values for the scattering 
length and the resulting xPT fit curve are shown in Fig. 2. In the same figure, 
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Fig. 2. Chiral extrapolation for the 1=2 pion-pion scattering length. The results in 
this work are shown together with the lattice calculations of NPLQCD [37, 38] and 
CP-PACS [36] and the direct measurement from NA48/2 at CERN [7]. 

we also provide a comparison to the lattice results of NPLQCD [37, 38] and 
CP-PACS [36] and the direct measurement from NA48/2 at CERN [7]. We 
find general agreement between our calculation and the results of NPLQCD 
at similar pion masses. In particular, the agreement between our results and 
NPLQCD suggests that the effect of the missing strange quark in our current 
calculation is small. Additionally, the ongoing effort of ETMC to include the 
dynamical effects of both the strange and charm quark [55, 56] will allow us 
to directly address this issue. 



To highlight the impact of the NLO terms in the xPT description of the pion 
mass dependence of m^a^^^ and to understand the role of yet higher order 
terms, we show the difference between the lattice calculations of the scattering 
length and the LO xPT prediction in Fig. 3. We find that the scattering 
lengths statistically agree with the LO xPT result for all lattice calculations 
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Fig. 3. Difference between tlie lattice calculation of the scattering lengths and the 
LO xPT prediction. The scattering lengths agree statistically with the LO xPT 
prediction for = 270 MeV to 485 MeV. 

with < 500 MeV. Accordingly, the NLO xPT functional form provides 
a reasonable description of the lattice results in the same region of m^r. As 
a further check, we fit our calculations to the NNLO form for mTrO^^^ [4, 54] 
and found m,rO^^^ = —0.041 (12) at the physical point. The statistical error 
is large, as one would expect given that our results already agree statistically 
with the LO xPT form, but the resulting NNLO extrapolation of m^ra^^^ does 
agree with the NLO fit. Given the size of the statistical errors, we are unable to 
make any meaningful estimate of the NNLO LECs, however, the effects from 
truncating the xPT series to NLO is included in our estimate of systematic 
errors. 



The systematic error on the extrapolated value of m^a^^^ and Z^^^ has sev- 
eral components. First, the systematic errors of the m^a^^^ that we obtain for 
each ensemble are propagated through the chiral extrapolation. This is accom- 
plished by again collecting all fit ranges for each ensemble with x^/dof < 2 
as earlier. This gives approximately 10^° xPT fits from which we randomly 
choose 2000 to sample the distribution of the extrapolated values of m^a^^^. 
As for the individual m^a^^^, we use the distribution method to determine 
an estimate of the systematic error due to the fit ranges from each ensemble. 
The second systematic uncertainty arises from the chiral fit itself. This is esti- 
mated by taking the difference in the extrapolated values from the NLO xPT 
fit to all six and just the lightest five ensembles. Finally, the extrapolation 
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to the physical point requires the experimental value for mTr//^. The exper- 
imental error on this quantity introduces an error that is nearly 50% of the 
corresponding statistical error and hence is also included. All three effects are 
added in quadrature to form the total estimated systematic error. Using the 
latest PDG [57] values of = 139.5702(4) MeV and /^+ = 130.4(2) MeV 
to determine the physical limit, we obtain the final result 

mXf = -0.04385 (28) (38) and /^f (// = ^,phy) = 4.65 (.85)(1.07) . 

This agrees with the previously mentioned results: the lattice calculation from 
NPLQCD [37,38], the so-called CGL analysis [3,4] and the E865 [2] and 
NA48/2 [7] measurements and represents agreement among the experimental 
and theoretical determinations of m^a^^^ at the 1% level. 

This accuracy of 1% must be understood as a combined theoretical effort from 
lattice QCD and chiral perturbation theory. The quantity m^a^^^, as well as 
a^^^ itself, vanishes in the chiral limit. This significantly constrains the chiral 
extrapolation of m^a^^^. In particular, rriT.a^^'^ is uniquely predicted in terms 
of m.„/ fj, at LO and depends only on one unknown constant, at NLO. 
This makes the chiral extrapolation of lattice results particularly accurate. 
Thus the 6% to 11% accurate results for the lattice calculation extrapolate to 
a 1% accurate determination of m,ra^^^. The final result differs from the LO 
xPT prediction by only a few percent, but to illustrate the power of combining 
lattice QCD and xPT, we note that the difference between mTrfl^^^ and the LO 
result is 0.00173 (47), which represents a 3.7cr shift that is due to the inclusion 
of the NLO effects as determined by matching directly to lattice QCD. 

Chiral perturbation theory plays a strong role in obtaining m^ra^^^ accurately, 
but it alone can not determine Only in combination with the lattice re- 
sults can we calculate /^^^ = 4.65(1.37). This calculation is just under 30% 
accurate, however, as we explain shortly, this easily exceeds the experimental 
determinations of this quantity. The experimental and phenomenological re- 
sults for m^ra^^^ at the physical point can be converted into a result for /^^^ at 
NLO. A simple analysis gives the following for 11="^: 3.0 ±3.1 (CGL), 0.0 ±10.3 
(E865 with xPT) and 3.0 ±2.8 (NA48/2 with xPT). This comparison demon- 
strates the particular advantage of lattice calculations arising from the ability 
to vary the underlying quark masses of QCD. 



5 Conclusion 

We have calculated the s-wave pion-pion scattering length in the isospin 1 = 2 
channel using the two-flavor maximally twisted mass lattice QCD configura- 
tions from ETMC. The pion masses ranged from 270 MeV to 485 MeV and 
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the lattice spacing was a = 0.086 fm. A second lattice spacing of a = 0.067 fm 
was used to demonstrate the absence of large lattice artifacts. This is only a 
single check, but when combined with the fact that the calculation is accu- 
rate to O(a^) due to the properties of maximally twisted mass fermions, it 
suggests that the lattice spacing dependence is mild. Furthermore, discretiza- 
tion errors vanish from the I = 2, I3 = ±2 channel at NLO, as shown 
by twisted mass xPT, hence we extrapolated our results for the scatter- 
ing length to the physical limit using continuum xPT at NLO. We investi- 
gated various systematic effects, and we found for the scattering length at the 
physical point m^a^^^ = —0.04385 (28) (38) and for the low energy constant 
/^^^(/i = /vr.phy) = 4.65 (.85)(1.07). These results are in good agreement with 
the previous lattice calculation from NPLQCD, the experimental determina- 
tions from E865 at BNL and from NA48/2 at CERN and the CGL analysis 
using various theoretical and experimental inputs. 
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